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AB3*RACT 

Computations  Indicate  that  when  utilised  In  an  idealised 
ramjet  at  a  flight  Mach  number  of  2.5  and  an  altitude  of 
37,500  feet,  a  constant  area  axial  flow  nuclear  heat  exchanger 
will  be  expected  to  produee  9*500  pounds  of  thrust  with  an  air 
inlet  Maeh  number  of  0.2,  a  weighted  radial  teriperature  aver- 
age of  2000°!*  and  a  frontal  area  of  7.7  square  feet.  Molyb- 
denum alloy  with  a  dlsillolde  coating  was  selected  as  the 
container  material.   .Sandwiched  between  the  molybdenum  was  a 
fissionable  alloy  of  highly  enriched  uranium  and  zirconium. 
From  a  nuclear  aspect  the  unref lected  heat  exchanger  becomes 
critical  at  a  diameter  of  lj.0  Inches  and  a  length  of  60  inches 
with  a  total  weight  of  8,052  pounds  which  includes  3,622 
pounds  of  90  per  cent  enriched  uranium. 
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That  a  constant-area  axial-flow  compact  nuclear  heat 
exchanger  as  envisioned  in  Figure  1  when  installed  in  an 
idealized  ramjet  would  produce  significant  thrust  was  the 
basis  of  this  thesis.   easibility  conclusions  were  limited 
to  the  information  revealed  by  heat  transfer  and  nuclear 
analysis  under  one  assumed  steady  state  condition*  The 
problems  of  nuclear  control,  shielding,  thermal  stresses, 
and  performance  under  different  inlet  Mach  numbers  were  not 
considered* 

The  nuclear  heat  exchanger  as  shown  in  figure  2  incor- 
porates thin  concentric  cylinders  which  are  the  solid  Tis- 
slonable  containers*   The  alloyed  fissionable  phase  is 
sandwiched  between  two  layers  of  non-fissionable  structural 
material  so  that  the  enriched  uranium  fuel  cores  are  sur- 
rounded by  a  continuous  network  of  structural  metal  for 
strength  and  containment  at  high  temperatures*  Air  flows 
within  the  annuli  and  acts  as  the  coolant  and  thrust  agent* 
It  is  recommended  that  the  nuclear  flux  be  flattened  by  a 
central  control  core  and  a  peripheral  reflector;  however,  the 
beneficial  effects  of  a  reflector  and  central  control  core 
were  not  considered* 

Prom  the  assumed  inlet  condition,  various  maximum  cyl- 
inder or  wall  temperatures  were  selected  arbitrarily  for 
thrust  results*   That  the  high  temperature  fuel  elements  and 
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structural  materials  were  possible  metallurgically  were  next 
investigated.  Whether  or  not  an  unreflected  homogeneous 
reactor  of  this  configuration  can  be  made  critical  was  deter- 
mined by  using  one  group  diffusion  theory. 


'inJv:->  mam  ?•*•% 


LIST  OP  3YKB0LS 


p 
A      Cross  section  area  normal  to  flow,  (ft  ),  or  the  mass 

number  of  an  element  (gm/gro  atomic  weight). 

A^     Wetted  area,  (ft2). 

•2 
I      Buckling ,  (cm  ). 

c  Local  velocity  of  sound,  (ft/sec). 

G  Specific  heat  at  constant  pressure,  (BTtJ/lb  *F). 

D  Diffusion  coefficient  for  flux,  (cm). 

D  Hydraulie  diameter,  D^  •  P,,  (ft). 

D,  Inside  diameter,  (ft). 

Dp  Outside  diameter,  (ft). 

-W  $  ffeutron  leakage,  (n/ctar   see). 

f  rlction  factor  defined  as  2c  T  J  pv2, 

g  Acceleration  of  gravity  (32.2  ft/sec  ). 

0  Mass  velocity  (lb/sec  f  t  )  • 

h  Coefficient  of  heat  transfer  between  fluid  and  surface, 

(BTU/hr  ft2  °P). 

he  based  on   (T^w  -  T0fMf>. 

hr  based  on  logarithm  swan  of   (T     .,  -  T) . 
u  ot  w 

h.    based  on  (T     M  -  T). 
■  otw 

h„  based  on  (T     „-  T  ) . 

3  OtW  O 

H  Extrapolated  distance,    (em). 

k  Neutron  multiplication  factor. 

K  Thermal  conductivity,    (BTtf/hr  ft2  «P/ft). 


. 


k 

L  Length  or  heat  exchanger  heating  element  In  the 

longitudinal  direction,  (ft). 
I     Haas  rate  of  flow,  (lb/see). 
M      Mach  number  or  molecular  weight  (graa/gm  molecular 

weight)* 
n      dumber  of  neutrons  per  car. 

V  ffumber  of  nuclei  per   cmr. 

23 
|L     Avogadros  number  (6.02  x  10  ^  nuclei/gm  atomic  weight). 

fLm  Prandtl  number ,  C  u/K. 

§mm         Reynolds  number,  ?v®e/V" 

!?«-    Stanton  number,  h Vc  l 
ox  e  p 

P  Absolute  stream  pressure  (lb/in  ) . 

r0  Isentropic  stagnation  pressure  (lb/in  ). 

^  Heat-transfer  rate  (BTU/hr). 

Q  Volumetric  heat  source  which  is  equal  to  FgJ, 

(BTU/seo  cur) . 

r  Radius  (ft). 

T  Absolute  stream  temperature  (*B). 

Th  Thrust  (lb). 

T*  Stagnation  temperature  (*R). 

TQ  Adlabatlc  stagnation  wall  teraperature  (*R). 
70  w   Stagnation  wall  temperature  (*B). 

▼  Velocity  (ft/sec). 

V  Volume  (ft^  or  cm?). 
M  Weight  (lb). 

x  Distance  along  heat  exchanger  heating  elements  (ft). 
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Subscripts 

b      Bulk  (T). 

f      I  11b  tf  =  (1^  t  tw)/2. 

w      Wall. 

0,1,2.     Station  references 

Greek 

a      Definition  of  vjl/L. 

t  Ratio  or  specific  heata  of  gaaf   Cp/Cv  a  1#^* 

ah     Change  in  enthalpy,     /lb). 

7p     number  of  neutron*  produced  per  absorption  in  the  fis- 
sionable isotope* 

X  Definition  of  W^f^3* 

A*    Transport  mean  free  path,  (cm), 

t*      Absolut©  viacosity  of  the  fluid,  (lb/sec  ft). 

]T     Average  cosine  of  the  scattering  angle  per  neutron 
collision  in  the  laboratory  aystem. 

JS  8maber  of  neutrons  produced  per  f iasion  in  the  fission- 

able isotope. 

p     Density,  (lb/ft3). 

2. 
<7~     Microscopic  croas  section  (cm  ) 

<T_  microscopic  absorption  cross  section  (era  ). 

a 

(T  microscopic  scattering  croas  section  (era  ). 
<T\  microscopic  fisaion  croas  section  (cm  ). 
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£  Macroscopic  croas  section  (on""  )f  (K<r). 

S  macroscopic  absorption  croas  section  (cm*  ). 
g  macroscopic  scattering  cross  seetion  (cm  )• 
£f   macroscopic  fission  cross  section  (cm*  ). 

S  0       tfumber  of  neutron  absorptions  /  oar  see* 

Tw     Tractive  force  at  wall  due  to  fluid  friotion  (lb/ft2). 

0  Flux  in  neutrons  /  cm  aec. 
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tiaplacian  operator. 
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LIMITING  CONDITIONS 

To  obtain  a  realistic  operating  condition  Tor  the  air- 
cooled  nuclear  heat  exchanger  the  following  conditions  were 
assumed: 


MQ  «  2.5,  h  «  37,500  ft,  T  ■  392#F,   P  «  3*09  lb/in2 


Godsey  and  Young  (1,  p. 221)  stated  that  the  total  temperature 
through  the  inlet  diffuaer  depends  only  on  flight  Each  num- 
ber and  ambient  temperature  and  Is  a  constant  regardless  of 
diffuser  losses  If  the  process  Is  considered  adiabatic. 
Therefore  the  relation 

T0.0  -  T0,l  •  T0.2  °  T(1  +  "^  M02)  (1) 

say  be  used  where  the  station  numbers  correspond  to  those  of 
ore  1.   Godsey  and  Young  (1,  p.221)  stated  that  the  total 
pressure  at  the  heat  exchanger  inlet,  I'   D,  was  dependent  upon 
the  flight  Kach  number,  ambient  pressure,  free  stream  diffu- 
sion losses,  duct  friction  losses,  and  duct  diffusion  losses. 
The  last  three  effects  may  be  grouped  under  the  term  ram  ef- 
fectiveness or  pressure  recovery  effectiveness,   ^rom  data 
presented  by  Godsey  and  Young  (1,  p»86),  the  ram  effectiveness 
at  Hq  »  2.5  was  given  as  approximately  0.98;  however,  subsonic 
diffusion  was  assumed  to  be  loss  free.  For  a  more  realistic 
condition  a  value  of  0.90  was  chosen.   The  above  ratios  may  be 


3 

written  out  In  the  following  expression  for  stagnation  pres« 
sure  at  the  heat  exchanger  inlet: 


P    .  .  fl  •  JL=l  m02)  V**  ,  V  actual         (la) 
0,2    0     ^   0         P0>2  ldeal 


Solutions  to  Equations  1  and  la  give  TQ  2  a  80l|°R  and 

P  -  ■  lj.6.5?  lb  per  sq  in.  for  inlet  conditions  to  the  heat 

exchanger • 

-lor  to  converting  the  above  inlet  stagnation  pressure 
to  stream  pressure,  it  is  necessary  to  examine  the  limiting 
conditions  for  heat  exchan  o.  Assuming  a  uniform  cross 
sectional  heat  exchanger  area  under  a  ^iven  inlet  stagnation 
temperature,  which  in  this  case  is  881f*R,  and  various  inlet 
Mach  numbers,  it  is  desired  to  exchange  heat  to  the  extent 
that  the  stagnation  temperature  at  the  exit  of  the  heat  ex- 
changer be  raised  to  a  maximum  value.  By  use  of  methods  out- 
lined in  Reference  (2,  p*l^#2U)  the  results  of  various  inlet 
Wach  numbers  are  shown  on  Figure  3  which  reveal  that  when  a 
certain  maximum  temperature  occurs,  thermal  choking  even 
under  frietionless  conditions  exists  in  the  heating  process 
and  definitely  limits  the  capacity  of  the  constant  area 
combustor  or  heat  exchanger.   3tated  another  way,  if  at  any 
position  downstream  for  a  constant  area  duct,  the  stagnation 
temperature  is  equal  to  T0  ^  maximum,  the  flow  is  choked 
since  the  Mach  number  is  unity.   Any  addition  of  friction 


will  reduce  the  downstream  Interval  for  TQ  -  maximum  to  occur. 
With  M2  m   0.2  and  TQ  ,  maximum  achieved,  $100*Ti$   the  stagna- 
tion  pressure  would  suffer  a  20  per  cent  loss  due  entirely  to 
the  heating  effect. 

Without  becoming  involved  with  a  specific  shock  cone  and 
diffuser  design,  the  inlet  Jfcch  number  was  set  at  0*2,  which 
resulted  in  the  following  stream  conditions: 


P2  -  ItS.l  lb/in2 


92  »  .1385  lb/ft3 
v2  *  291  ft/sec 
T2  «  880*K 
G  *  \\Q*h   lb/ft   sec 


As  shown  by  Pigure  2,  the  heat  exchanger  elements  are 
thin  hollow  constant  concentric  area  right  cylinders  which 
are  adaptable  to  one-dimensional  compressible  gas  flow  treat- 
aent  without  introducing  serious  errors.   This  assumes  that 
in  any  one  annulus  the  change  of  stream  properties  in  the 
direction  of  flow  is  much  larger  than  in  the  radial  direction. 


.    . 
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HEAT  TRANSFER  AMTX&13 

WcAdama  (3,  p. 289)  correlated  the  experimental  data  of 
Kays  and  London  (k)   and  baoauae  compactness  of  an  air-borne 
heat  exchanger  ia  an  important  parameter,  a  comparison  of 
heat  tranafer  coefficient a  per  unit  volume  for  a  given  fric- 
tion power  per  unit  volume  showed  ruffled  fine  and  louvered 
plate  fine  as  being  the  moat  effective.   It  appears  that  the 
highest  ratio  of  heat  transfer  to  pumping  power  ia  obtained 
when  form  drag  la  absent,  or  when  all  the  dra£  ia  due  to  akin 
friction  which  ia  the  caae  in  common  for  flat  platea  or  fins* 
Even  though  externally  finned  tubes  are  of  value  in  increas- 
ing the  rate  of  heat  tranafer,  friction  loasea  become  pro- 
hibitive for  this  type  of  inatallatlon.  i'or  these  reasons 
thin  concentric  hollow  cylinders  were  chosen  as  the  fuel  ele- 
ments  of  the  heat  exchanger  and  were  restricted  to  constant 
longitudinal  crosa  section  area  for  nuclear  simplification. 

Development 

If  heat  loss  and  longitudinal  thermal  conductivity  are 
neglected,  the  rise  in  the  enthalpy  of  the  air  per  second  ia 
related  to  the  effective  coefficient  of  heat  transfer  from 
an  annulua  by  a  differential  equation  which  is  developed  in 
Appendix  I  as 


ii#  ' 
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I   <D22  -  D^lpT  Cp  <JTQ  -  he«(D2  ♦  D1)(t0#v  -  TQ#av)    «x.      (2) 

The  particular  temperature  differential  is  shown  schematically 

in  Figure  lj..      The  temperature   increment  for  heat  transfer  is 

the  difference  between  the  actual  wall  temperature,   l-     ,   and 

»w 

the  temperature  Tn       ,   which  the  wall  would  assume  for  the 

u,  aw 

same  flow  conditions  but  with  zero  heat  flux*   Mc Adams  et  al« 

(5)  found  that  the  effective  heat  transfer  coefficient,  h  , 

expressed  in  terms  of  the  Stanton  number,  N,,-,  was  independent 

of  the  temperature  potential,  TQ   -  TQ,  when  based  on 

*n  «  •  ^n  «,,  while  the  coefficients  h  and  h  were  not  inde- 
ed   o,aw  ■      s 

pendent. 

The  Reynolds  analogy  with  NpR  «  1  relates  the   Stanton 
number  to   the  friction  factor  by 

V       f 

•pr-g  ■  *g  •  (3) 

P 

Prom  experimental  data  for  a  moderate   AT,    h/n     »  60,  and 

e 

S,000<MRS<  200,000,   the  following  equations  from  MoAdams 
(3,   pp •155,219)   may  be   combined 


•2 


f  »   .<W(D0/|ifr*  (k) 


A  (¥))      -  T^TT?  <*> 


pbw     \  *  /f  (D0/|tf ) 


to  produce 

w  ''''"Tf  •  — -w—  (6) 

where,  by  definition,  Az   (Kp^)?         •  Since  Equations  3  and  6 
are  essentially  the  same,  the  limitation  to  a  moderate  tempera- 
ture change  may  be  removed  without  introducing  serious  errora 
when  using  the  previously  defined  temperature  potential.   The 
following  relationship,  where  specific  heat  is  considered  a 
constant,  Is  used  for  the  remainder  of  the  analysis: 


h 


r-§-"-r  (7) 


The  friction  factor  of  Equation  7  is  set  by  Equation  Ij.  which 
Is  dependent  upon  Re;/nolds  number  and  the  degree  of  roughness. 
Since  mass  velocity,  G,  and  specific  heat  are  essentially 
constant,  the  friction  factor  may  be  considered  to  depend  upon 
hydraulic  diameter,  viscosity,  and  degree  of  roughness*  An 
average  value  of  /If  »  *006  was  obtained  from  Equation  I4.  for 
1/2* <  De<l-1/2W  and  ll*00°R  <  p.f  <1600°H.   The  variance  of  the 
friction  factor  within  the  above  range  of  hydraulic  diameter 
may  be  compensated  for  by  added  degrees  of  roughness;  however, 
Cope  (6)  found  that  with  water  in  the  turbulent  region  even 
though  friction  was  six  times  that  for  smooth  pipes,  the  heat 
transfer  was  only  20  to  100  per  cent  greater  than  for  smooth 
tubes*  MeAdams  (3,  p»223)   states  that  because  of  a  very 
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small  temperature  change  the  heat  transfer  coeff icients  ob- 
tained were  somewhat  uncertain.  Correlation  with  air  may  not 
apply  due  to  different  Prandtl  number  effect j  however,  there 
is  a  possibility  that  with  small  variation  in  hydraulic  diam- 
eter the  friction  factor  could  remain  a  constant  by  varying 
the  degrees  of  roughness, 

Experimental  data  have  shown  that  for  turbulent  flow  in 
rectangular  ducts  and  annular  spaces  friction  may  be  evaluated 
from  data  for  circular  pipes  using  a  hydraulic  diameter  shown 
in  Figure  2  as 


De  *  ^h  *  ~^(b2  ♦  D*>   ■  D2  -  Dl   •    <8> 


A  combination  of  Equations  2  and  7  with  Equation  8  results  in 

4j!o-  I&s-Is**  aargs  .   (9) 

If  a  gas  stream  of  uniform  temperature  is  brought  to  rest 
isentropically,  as  at  the  stagnation  of  a  blunt  body,  the  tem- 
perature rise  for  an  ideal  ,*ras  will  be 


T(l  «  J^=l  M2)      .    (10) 


For  convenience  a  temperature  recovery  factor  N^  as  shown  in 
Figure  1*  is  defined  as 


11* 


I     -  T    T     •  T 

Of** o^aif .     (a) 

T0  -  T    V  /2%   JC 

Mr 


By  a  combination  of  Equation  10  and  11,  the  stream  temperature, 
T,  may  be  eliminated  to  produce 


'g.a»  ,  ,   <*-%>  4^**  ,12) 


Substitution  of  Equation  12  into  9,  yielda 


dT 


0 


f^r"1* — ,  A  * -1  „g —     2Afr     •    (13) 


T0        0  1  ♦  4*=*  H2 


r 

The  generalised  differential  Maoh  number  equation  of 
Reference   (2#  P.I4..I2)   aimplifiea  to 


(1  ♦    JT^Ml  ♦  -J|=i  M2)  dTQ 

*0 


M2*!  *  JLd.  M2)  ., 

S kf  S*  (1U) 

1  -  ?r  ue 


with  negligible  structure  drag  forces  and  constant  specific 
heat,  molecular  weight,  area  and  mass  flow  rate. 


' 
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Equation  13  was  substituted  Into  Equation  Ik  which  gave 


A* m    rr? — 


I 


(1  +  *  M2)/* 


"O.w 


-  1  ♦ 


l-H^^^H2 


1  ♦ 


+  2r 


(15) 


For  a  turbulent  boundary  layer  over  a  flat  plate, 
MoAdams  (3,  p. 311)  recommends  that  the  temperature  recovery 
factor  and  the  Prandtl  number  be  related  by 


W      a,  M      1/3 

aRF    "PR 


(16) 


MeAdams  et  al.  (5)  reported  recovery  factors  for  subsonic  flow 
of  air  near  the  outlet  and  of  a  long  smooth  tube  were  in  good 
agreement  with  those  predicted  by  Equation  16.  Appendix  A 
shows  the  procedure  for  solving  Equations  13  and  15  for  one 
Mach  number  and  stagnation  temperature  at  an  assumed  position 
downstream*  From  these  calculations  it  was  evident  that  the 
contribution  of  the  recovery  factor  term  &**&»   wa8  negligible 
for  Mach  numbers  up  to  approximately  0*8}  however,  the  re- 
sulting stagnation  temperatures  will  be  nearly  1.5  V*r   cent 
higher  than  those  which  have  Incorporated  this  term. 

W   the  remainder  of  this  analysis  the  recovery  factor 
and  Prandtl  number  were  assumed  unity  which  reduced  Equation 
13  to 


16 


■ 

*0.w   . 

T   " 


Zt 


&  (17) 


It  is  of  interest  to  not©  that  Equation  17  oould  bo  expected 
from  experimental  data  by  MoAdams  (3,  p.31ij.)  because  the  heat 
transfer  coefficient,  h  ,  decreases  and  approaches  the  constant 

value  of  h  for  temperature  differences  greater  than  100°P. 

e 

For  a  high  temperature  potential,  h  approaches  h  and  may  be 
used  directly  in  Equation  1. 

Equations  II4.  and  17  were  integrated  approximately  over 
ahort  intervals  of  Kach  numbers  by  using  the  procedure  out- 
lined in  Reference  (2,  I4..6O).  Appendix  B  contains  the  step- 
wise procedure  used  and  Tables  1-3,  Appendix  B,  have  the  com- 
plete data  for  stagnation  and  stream  properties  as  a  function 

of  L/D  ratios  for  various  longitudinally  constant  wall  tew 

e 

peratures.   Stagnation  temperatures  and  pressures  versus  L/D 
ratios  are  shown  in  Figure  £• 

wever,  from  a  nuclear  standpoint  the  longitudinal  flux 
of  fuel  element  cylinders  will  follow  a  sine  distribution  with 
the  maximum  flux  produced  at  the  midpoint.   This  does  not 
imply  that  the  maximum  wall  temperature  will  occur  at  the 
longitudinal  midpoint.  By  use  of 

t    -  t    1    O.w  max  -  0.2*     MA* 
*      '      1  -  sec  a  max 


r 


•- 
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which  was  developed  In  Appendix  C  end  by  using  the  previously 
•elected  temperatures  for  maximm  wall  temperatures,  the 
stagnation  termers tures  versus  V^«  ratios  were  computed  as 
shown  in  Table  k   of  Appendix  0  and  plotted  in  lure  5»   Stag- 
nation pressures  for  variable  wall  temperatures  were  obtained 
by  interpolation  from  previously  calculated  pressures* 

An  isentropic  expansion  using  air  data  from  Keenan  and 
Kaye  (7)  was  assumed  in  order  to  obtain  the  change  of  enthalpy 
from  stage  3  to  k  of  Figure  !•   The  Brayton  cycle  as  shown  in 

Figure  6  represents  an  overall  view;  however,  at  various  Vb 

e 

ratios,  the  chanre  of  enthalpy  was  computed  using  an  ambient 
pressure  of   3*09  lb  per  so  in*  Figure  7  shows  the  variance 
of  enthalpy  for  both  constant  and  variable  wall  temperatures* 
The  right  hand  ordinate  of  ?igurs  7  related  change  of  enthalpy 
to  pounds  of  thrust  per  sq  ft  by 


Th  ..  0 


-r m  f  {\m  V      ■  (19) 


where  vQ  *  2.5c  and  ▼•    m  223.7Vn~0  ^  -  ii. 


Discussion 


The  calculated  values  of  stagnation  temperatures  and 
pressure  as  shown  in  Figure  5  do  not  show  the  contributing 
effects  on  heat  transfer  and  friction  caused  by  the  length 
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to  diameter  ratio  and  the  effects  of  sudden  enlargement  or 
contraction  of  flow.   The  form»T   effect  is  small  since  the 
L/D  ratio  for  this  heat  exchanger  should  be  greater  than  60; 
the  latter  effect  may  still  be  appreciable  even  though  the 
peripheral  thickness  of  the  cylinders  are  in  the  order  of  0.0< 
to  0.07  in.   The  specific  heat  ratio  was  set  at  l*k»   and  the 
friction  factor,  specific  heat,  and  mass  velocity  were  con- 
stant throughout  the  heating  process.  The  cosine  flux  distri- 
bution in  the  longitudinal  direction  is  justifiable  for  a 
finite  cylinder. 

The  temperatures  obtained  by  normal  forced  convection 
Methods  for  an  assumed  TQ   ■  21^6 0°R  are  shown  in  Figure  5>  for 
comparison*   The  temperatures  show  good  correlation  with  stag- 
nation temperatures  obtained  by  the  Reynolds  analogy  but  are 
considerably  higher  than  the  stream  temperatures.   The  forced 
convection  formulas  will  not  show  any  decreasing  trend  in 
stream  temperatures  or  a  maximum  stagnation  temperature  as 
heat  is  continually  added  as  will  the  Reynolds  method  for  com- 
pressible fluid  which  will  eventually  produce  choked  flow  at 

a  certain  r/D  ratio  downstream, 
e 

Figure  7  shows  that  the  effect  on  thrust  reduction  was 
quite  teyere   for  a  2000*F  constant  wall  temperature  versus  a 
2000°F  maximum  wall  te?r$>erature;  the  thrust  per  square  foot 
was  reduced  by  18  per   cent.  For  a  uniform  velocity  at  the 
exit  of  stage  U,  a  constant  change  in  enthalpy  is  required. 
If  the  neutron  flux  was  such  as  to  produce  a  series  of  maximum 
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wall  temperatures  from  1600  to  220  °P  and  with  a  constant 
enthalpy  change  of  200  BTU  per  lb  imposed,  the  L/D^  ratio 
would  vary  t9  t   to  100.   MJfcM  all  cylinders  are  desired 

to  have  the  same  length  for  nuclear  considerations,  the  hy- 
draulic diameters  would  vary  by  fc?  per  cent  with  decreasing 
hydraulie  diameters  ss  radial  distance  is  increase  .     is 
diameter  variance  would  still  be  subject  to  the  limitation 
imposed  on  the  hydraulic  diameter  due  to  a  constant  friction 
factor. 

To  increase  the  thrust  which  is  based  on  the  heat  ex- 
changer ei'os  a -sectional  area,  each  fuel  cylinder  should  oper- 
ate at  the  maximum  enthalpy  change  as  shown  by  Figure  7  and 
the  thrust  would  be  a  weighted  average  depending  on  the 
frontal  area  of  the  heat  exchanger.   With  a  constant  length 
assumed  the  hydraulic  diameter  would  vary  by  only  12  jser   cent. 
However,  each  annulua  would  be  operating  at  a  different  veloc- 
ity level  which  may  create  internal  vibration  within  the  ram- 
jet pod.   With  a  weighted  radial  average  of  2000°F  for  the 
cylinder  temperatures  encompassing  7.7  sq  ft  of  heat  exchanger 
cross-sectional  area,  from  Figure  7  the  resultant  thrust  would 
be  approximately  9,500  lb  at  a  flight  Ma eh  number  of  Mmt   at 
37,500  ft. 

With  inlet  wash  numbers  greater  than  0.2,  the  E/P  ratio 

e 

would  be  substantially  reduced.      i*rom  Reference    (2,   p»li31), 
with  K2  «  0.5,    f  ■  0.005,    TQ  2  ■  10O0«R,    P0  2  -  10  lb/in2,   and 
T0  w  *  2000«B  &***  1500*R,    choked   flow  occurred  at  an  L/D     ratio 
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of  20  and  35,  respectively.  To  resist  the  pancaking  effect 
which  »iay  reach  the  limit  for  nuclear  criticality,  the  inlet 
Mach  number  appears  to  have  a  definite  upper  limit  when  heal 
lug  is  desired. 

Conclusion 

For  the  assumed  inlet  conditions,  cylindrical  fuel  ele- 
aents  and  wall  temperatures ,  the  amount  of  heat  transferred 
was  substantial  enough  to  produce  1,000  to  1,1^00  pounds  of 
thrust  per  sq  ft  if  heat  exchanger  frontal  area  with  a  maximum 
wall  temperatures  ran^-e  from  1600  to  2200#t*  at  a  flight  Mach 
number  of  2.5  and  an  altitude  of  37,500  ft.  With  a  weighted 
radial  average  of  2000°F  for  the  cylinder  temperatures  encom- 
passing 7.7  sq  ft  of  heat  exchanger  area  the  resultant  thrust 
was  9,500  lb. 


if 
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Ml  COSSIDFRATIOH 

Due  to  the  absence  of  effective  moderation  within  this 
nuclear  heat  exchanger,  the  unit  may  be  considered  as  a  fast 
neutron  reactor  and  aa  euch  the  restrictions  normally  imposed 
on  structural  materials  by  high  thermal  neutron  absorption 
may  be  removed.   Whether  any  structural  material  now  in  exiat- 
•nce  can  be  used  to  contain,  protect,  and  support  the  fission- 
able phase  for  the  high  temperatures  required  is  the  subject 
of  the  following  development  and  discussion. 

Development 

The  nuclear  heat  exchanger  as  shown  in  Figure  2  incor- 
porates thin  concentric  cylinders  which  are  the  solid  fission- 
able containers*   The  enriched  uranium  fuel  cores  are  to  re- 
main a  solid  to  prevent  serious  containment  problems.  The 
thin  sheets  of  the  fissionable  phase  are  sandwiched  between 
the  inner  and  outer  layers  of  the  non-fissionable  structural 
metal  which  predominates  in  volume  so  that  the  alloyed  en- 
riched uranium  fuel  sheets  are  surrounded  by  a  continuous 
network  of  structural  metal  for  strength  at  high  temperatures, 

The  following  are  some  criteria  that  were  considered  for 
the  selection  of  fuel  elements  and  structural  material: 
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A.  Fuel  element  criteria 

1.   Melting  point  above  2500°P 
neutron's  energy  maintained 

3,  Mechanical  properties  similar  to  the  container 
material  during  fabrication  and  high  tempera- 
ture operation. 

4.  High  thermal  conductivity  and  low  specific 
heat  for  good  heat  transfer 

£.   l-ase  of  fabrication 
.   Container  material  criteria 

1#   High  temperature  mechanical  properties  under 

static  and  small  radial  dynamic  loads 
2*  Resistance  to  oxidation  at  high  temperatures 
3*   High  resistance  to  thermal  shock 
\*      Neutron  energy  maintained 

Ease  of  fabrication,  rolling,  and  weldabllity. 

riscuselon 

The  melting  point  of  natural  uranium  is  2072*^;  for  the 
required  temperature  range  of  2CK  0  to  2200* ?$   the  unalloyed 
uranium  would  become  a  liquid  and  present  serious  containment 
problems*   In  order  to  increase  the  melting  point,  the  various 
binary  diagrams  of  uranium  alloys  of  Sailer  (8,  p,ijl3-lAt*-) 
revealed  that  a  uranium-zirconium  system  of  20.7  pT   cent  zir- 
conium by  atomic  weight  would  have  a  melting  point  of  Zk$0*F • 
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Elements  with  lees  than  a  mass  number  of  kO   were  not  consid- 
ered because  of  the  Increased  slowing  down  effect  on  the 
neutrons. 

Zirconium  alloys  at  high  temperatures  suffer  a  loss  of 
strength  and  a  loss  of  corrosion  resistance.  As  reported  by 
Dayton  (9,  p.l*60) ,  the  strongest  alloy  tested  at  930*F  creeps 
at  a  rate  of  10~^  per  cent  per  hour  with  a  stress  of  only 
I4.OOO  lb  per  sq  in.   Since  a  thin  sheet  is  sandwiched  between 
two  Isyers  of  structural  material,  the  loss  of  strength  and 
corrosion  resistance  are  secondary  to  the  requirement  that  the 
fissionable  material  remain  a  solid*   The  coefficient  of 
linear  thermal  expansion  for  20.7  per  cent  zirconium  by  weight 
was  not  available;  however,  Dayton  (9.  p.  1*6)4.)  reported  a  co- 
efficient for  pure  zirconium  at  932*P  for  the  "a"  axis  of 
5*8  x  10   while  that  of  pure  uranium  was  anisotropic  as  shown 
by  Sailer  (8,  p«389)«  The  coefficient  of  thermal  expansion 
for  this  particular  alloy  was  not  found;  however,  from  Sailer 
(8,  pm23k)   the  composition  for  20#7  per  cent  zirconium  by 
weight  indicated  3  phase  changes  for  the  temperature  range  up 
to  2200°P  and  whether  the  anisotropic  effect  of  uranium  would 
result  in  volume  changes  so  severe   as  to  make  bonding  or 
cladding  impractical  was  not  determined* 

The  low  thermal  conductivities  of  both  zirconium  and 
uranium  are  partially  cowpensated  for  by  low  specific  heats 
and  the  thin  sheet  effect  of  high  area  to  volume  ratio. 
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Uranium  can  be  hot  rolled  without  danger  of  cracking. 
Whether  this  particular  alloy  can  be  rolled  Into  sheets  with 
thicknesses  as  small  as  0*03  inch  was  not  determined. 

The  container  material  for  temperatures  up  to  2200 *F 
presents  a  considerable  problem  when  restricted  to  present  day 
material  availability.  The  wrought  cobalt-base  alloys,  such 
as  Haynes  3tellitc-25>,  Allegheny- Ludlum  Steel  Corporation 
3-£9Q  and  3-816,  and  Westinghouse  Electric  Corporation  Re- 
fractaloy  26  and  70,  appear  to  have  excellent  corrosion  re- 
sistance up  to  temperatures  as  high  as  2000*?,  Eichen  and 
Jackson  (10,  p«55>6) ,  in  an  oxidizing  atmosphere*   However, 
Riehen  and  Jackson  (10,  p»j>39)  reported  the  ultimate  strength 
of  Haynes  Alloy  25   at  2200*F  to  be  lj.,000  lb  per  sq  in  for  a 
sheet  0.0U2  inch  in  thickness  and  annealed  at  2200°:  .   Another 
major  drawback  may  be  the  high  residual  activity  in  the  form 
of  gamma  radiation  which  is  characteristic  of  cobalt-base  al- 
loys. 

Molybdenum  is  of  interest  as  a  high-temperature  material 
because  of  its  high  melting  point,  1^750*F,  high  strength  and 
stiffness  at  elevated  temperatures,   tforthcott  (11,  p,106) 
reported  creep-rupture  strength  for  wrought  arc  cast  molyb- 
denum, .U5  per  cent  titanium,  of  25,000  lb  per  sq  in  at  2000*? 
for  100  hours.  The  thermal  expansion  coefficient  as  reported 
by  Horthcott  (11,  p. 26)  ranged  from  5»1  x  10~*  at  932°F  to 
7.2  x  10   at  3632* '.  Thus  the  thermal  expansion  of  both  pure 
sirconium  and  molybdenum  alloys  appear  to  be  very  similar. 
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However,  alloying  with  uranium  might  alter  the  thermal  ex- 
pansion characteristic  considerably. 

Uneoated  molybdenum  is  useless  as  a  Material  of  construc- 
tion in  oxidizing  atmospheres  above  !060*P«   Norden  (12f 
pp.293-30lj.)  reported  that  so  far  no  molybdenum-base  alloys 
have  been  developed  with  satisfactory  oxidation  resistance; 
however,  Kiester  et  al.  (13)  reviewed  progress  lasde  with  al- 
loying techniques.  Sherwood  (114.,  p.26)  stated  that  molybdenusi 
specimens  protected  by   coatings  0*003  in  thick  and  consisting 
primarily  of  molybdenum  dlslliolde  have  successfully  resisted 
oxidation  at  3092*F  for  more  than  1000  hours.   However,  the 
coating  must  be  pore-free  in  order  to  prevent  serious  loss  of 
material.   Effects  of  the  disiliclde  coating  on  heat  transfer 
and  thermal  cycling  were  not  included. 

The  high  thermal  conductivity  of  molybdenum,  .32  cal  per 
sec  cm  *C,  and  low  speeific  heat,  .061  cal  per  gm   *Cf  permit 
it  to  be  heated  or  cooled  rapidly  with  low  resultant  stress 
which  give  it  excellent  resistance  to  thermal  shock.  The  heat 
transfer  characteristics  are  excellent  due  to  the  high  thermal 
conductivity.   The  neutron  moderation  effects  are  small  since 
it  has  a  mass  number  of  96. 

As  brought  out  by  Iforden  (12,  p*298)  a  chief  advantage 
of  molybdenum  over  ceramics,  cermets,  and  some  super  alloys  is 
the  variety  of  forms  and  sizes  available  on  a  production  basis. 
By  the  arc-cast  methods,  sheets  from  0.005  to  0.020  in.  thick 
and  36  in.  wide  are  available.  Molybdenum  has  been  joined 
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successfully  by  spot  welding,  arc  welding  soldering,  and 
brazing  but  precise  controls  must  be  used  since  excess  heat- 
ing will  result  in  extremely  brittle  welds. 

Conclusion 

Molybdenum  alloy,  0.1$  per  cent  titanium,  with  a  protec- 
tive disilielde  coating  appeared  to  have  the  best  probability 
of  withstanding  temperatures  up  to  22G0°F  in  an  oxidising  at- 
mosphere besides  being  able  to  withstand  thermal  shock  and 
offer  near  optimum  heat-transfer  characteristics.   'The  inner 
fissionable  phase  of  enriched  uranium  alloyed  with  20»7  per 
cent  by  weight  of  zirconium,  even  though  soft  and  pliable  at 
this  temperature  range,  which  could  be  an  advantage  as  far  as 
the  container  material  Is  concerned,  should  expel  the  heat 
produced  internally  due  to  the  high  surface  to  volume  ratio. 
Whether  or   not  the  thermal  expansion  coefficient  of  the 
uranium- zirconium  would  be  compatible  with  molybdenum  and  that 
a  metallurgical  bond  would  be  formed  with  molybdenum  was  not 
ascertained. 
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The  actual  heterogeneous  reactor  configuration  aa  shown 
in  Figure  2  was  treated  as  a  homogeneous  unref looted  system 
using  one-group  diffusion  theory  to  obtain  the  critical  mass 
of  uranium.  The  flattening  of  flux  and  redaction  of  critical 
mass  by  use  of  a  reflector  were  not  considered.  The  neutron 
flux  necessary  to  produoe  a  maximum  fuel  element  temperature 
was  calculated.  The  control  mechanism  necessary  to  maintain 
the  flux  level  was  not  considered. 

Development 

The  diffusion  equation  for  thermal  neutrons  for  steady 
state  conditions  as  developed  by  Glasstone  and  Edlund  (15, 
p. 192)  is 


DV2J*  (r)  -  fia  ft   (r)+  pq  (r,  T  th)  -  0  (J 


For  an  idealized  fast  reactor  composed  of  atoms  of  material 
so  heavy  that  neutrons  upon  colliding  rebound  without  losing 
an  appreciable  amount  of  energy,  the  slowing  down  density  term, 
pq  (r, Tth),  may  be  considered  negligible.   If  fast  neutrons 
occur  only  at  one  energy,  the  source  term  can  be  written  as 
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Using  Glasstone  and  Edlund*s  (15)  treatment  the  multi- 
plication factor  may  be  simplified  to 


Ke*  v 


by  assuming  that  the  fast  fission  factor,  the  resonance  escape 
probability,  and  the  fast  neutron  utilization  are  essentially 
unity*  The  fast  neutron  utilization  term  implies  that  the 
eapture  of  neutrons  by  material  other  than  the  fissionable 
phase  is  negligible*   Shapiro  (16,  p»130)  states  that  the  ra- 
diative capture  cross  sections  for  elements  having  mass  numbers 
less  than  100  and  neutrons  numbers  less  than  60  are  not  greater 
than  £0  millibarns*  If  the  fast  neutron  utilisation  factor 
were  not  unity,  the  treatment  eiven  by  Kaplan  (17,  p»525)  *or 
thermal  neutrons  would  still  apply. 

A  combination  of  Equations  20  to  22  yields 

V2/f  (F)  +  B2  0   (F)  -  0  (23) 


where  B2  i  Sft  (  ^  -  l)/D.  (21;) 


For  a  weak  absorber  as  described  by  Glasstone  and  Fdlund 
(15,  pp. 92, 98),  the  diffusion  coefficient  is  related  to  the 
transport  mean  free  path  by 
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p  m   -^p  *  4  t  (25) 


where   At  3  As  (1  -  f0)  and  £  •  2/3A.  (26) 


As  a  first  approximation  the  macroscopic  transport  and 
scattering  cross-sections  are  considered  to  be  essentially 
equal  except  in  cases  where  the  actual  values  are  known.   The 
atomic  masses  of  the  uranium,  zirconium,  and  molybdenum  are 
high  enough  to  make  the  |t0  term  of  Equation  26  negligible.  As 
a  second  approximation  the  elastic  scattering  cross  section  is 
expected  to  be  equal  roughly  to  the  total  cross  section  in  the 
middle  energy  region  of  1  to  $00  kev  where  neutron  gamma  cross 
sections  are  sraall,  Shapiro  (16,  p«78).  The  constants  used  in 
this  analysis  are  listed  in  Appendix  ■  with  additional  refer- 
ences* 

Substitution  of  Equation  2$  into  2k  results  in 

B2  m  3£tSa(7]  -  1)  (27) 

The  solution  to  Equation  23  for  the  buckling  in  terms  of 
spatial  dimensions  of  a  right  circular  cylinder  as  developed 
by  Olasstone  and  Edlund  (15#  p«213)  is 

B2  -  (^22)  '  ♦  (!)  2  (23) 

where  P  and  H  are  the  extrapolated  radius  and  height  of  the 
cylinder  in  centimeters. 


_A 


wUx&vqGi 


- 


•d: 


■ 


• 


30 
A  combination  of  Equations  27  and  28  yields 

3fit2a(7,.l)  -  (HF02+  (i)'  (29) 

Equation  29  is  the  approximate  critical  equation  for  an 
Idealized  fast  homogeneous  unref lectcd  reactor  in  the  form  of 
a  right  circular  cylinder.   Appendix  P  shows  the  calculations 
that  were  used  to  determine  the  critical  mass  of  uranium  for 
a  right  circular  cylinder  with  a  physical  diameter  and  length 
of  \\2   and  60  In,  respectively. 

The  various  component  weights  which  were  based  on  volume 
fractions  with  corresponding  densities  are  calculated  in  Ap- 
pendix D. 

For  the  flux  required  to  produce  the  maximum  temperature 
change,  the  relationship 

4        -,  (T0.3  '  T0.2>  °  Co  » 

'wax       UE£rt  (L/De)  (30) 

was  used.   The  result  is  shown  in  Appendix  D« 

Discussion 

The  length  of  the  heat  exchanger  was  limited  to  60  in. 
because  Figure  7  showed  that  the  maximum  thrust  per  sq  ft  was 
obtained  with  an  approximate  l/D^  ratio  of  120.  This  produced 
a  0.50  in.  hydraulic  diameter  or  a  0#25  in.  fuel  element  spaolng 
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The  uranium-zirconium  alloy  with  a  composition  of  79.3 
per  cent  enriched  uranium  by  weight  was  selected  because  of 
its  2 , \&0*F   melting  point.  A  larger  percentage  of  uraniua 
would  reduce  this  melting  point  which  is  considered  a  minimum 
for  the  high  temperatures  desired*   The  alloy  composition  was 
considered  a  constant  and  treated  as  a  mixture  for  volume 
fractions  since  the  actual  density  was  unknown* 

The  component  volume  fractions  of  the  heat  exchanger 
resulted  from  trial  and  error  thickness  selections  for  the 
uranium-zirconium  and  molybdenum  alloys  which  satisfied  Equa- 
tion 29  for  a  physical  diameter  and  length  equal  to  or  less 
than  kZ   and  60  in,  respectively*   Since  the  physioal  length 
was  fixed  previously  for  maximum  thrust  per  sq  ft,  the  radius 
was  considered  variable* 

The  nuclear  reactor,  when  treated  as  a  homogeneous  un- 
reflected  reactor  according  to  Equation  29,  was  determined  to 
be  critical  at  a  physical  diameter  of  !|0  in  and  a  length  of 
60  in.  The  large  critical  volume,  i^»6  cu  ft,  was  due  to  the 
low  macroscopic  absorption  of  uranium,  0*0056  cm  ,  which  re- 
sulted indirectly  from  the  75  per  cent  volume  of  air*   A  de- 
crease in  the  volume  fraction  of  uranium  would  decrease  both 
the  total  macroscopic  absorption  and  transport  cross  sections 
with  a  resultant  increase  in  radius  if  the  length  remained 
constant. 

The  weight  of  enriched,  90  per  cent  uranium  required  for 
the  critical  volume  of  U3*6  ou  ft  was  3,622  lbs.   The  amount 
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of  zirconium  and  molybdenum  was  960  and  3,ij7o  lbs,  respective- 
ly, which  resulted  In  a  total  heat  exchanger  weight  of  8,052 
lbs.  At  current  prices  of  approximately  30  dollars  per  gram, 
50  million  dollars  would  be  required  for  the  90  per  cent  en- 
riched uranium* 

The  flux  that  was  necessary  to  generate  a  maximum  wall 
temperature  of  2200*F  with  a  0*0315  inch  fuel  element  thick- 
ness was  8  x  10  ^  n  per  sq  cm  sec  which  assumed  sero  radial 
flux  depression  across  the  fuel  element. 

Since  the  radial  flux  will  follow  an  approximate  cosine 
distribution*  it  Is  apparent  that  for  an  average  maximum  wall 
temperature  of  2000°P  which  was  assumed  to  obtain  9,500  lbs 
of  thrust  with  a  7.1   sq  ft  frontal  area,  the  fuel  elements 
near  the  centerline  would  operate  at  temperatures  exceeding 
2200*P.  A  proposed  central  control  rod  and  reflector  would 
tend  to  flatten  the  flux;  however,  the  flux  envelope  resulting 
from  these  effects  was  not  investigated. 

The  application  of  the  one-energy  group  diffusion  theory 
to  this  analysis  was  a  realistic  approximation  because  the 
transport  mean  free  path,  8»7  cm,  was  small  in  comparison  to 
the  boundary  distance.   The  extrapolation  distance,  R-r,  is 
not  generally  equal  to  0.71  A.  for  a  fast  reactor  and  its 
exact  value  can  be  obtained  only  by  complicated  methods,  Glas- 
stone  {20,   p.217).   However,  as  pointed  out  by  Glasstone  (20, 
p. 217),  there  is  probably  no  point  in  a  faat  reactor  that  is 
more  than  two  scattering  or  transport  mean  free  paths  from  a 
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boundary.  This  statement  is  supported  by  the  large  quantities 
of  fuel  required  for  the  heat  exchanger  which  has  a  minimum  of 
11  transport  mean  free  paths  in  the  radial  direction. 

An  Intermediate  neutron  energy  reactor  may  offer  the  pos- 
sibility of  reducing  the  fuel  quantity  and  weight  while  still 
Maintaining  a  large  volume  fraction  of  air  for  compact  heat 
transfer.  The  slowing  down  density  term  of  Equation  20  was 
neglected  because  air,  which  composed  75  per  cent  of  the  heat 
exchanger  volume,  has  a  negligible  moderation  effect  on  neu- 
trons. To  continue  this  simplification  the  material  for  the 
fuel  alloy  and  container  was  selected  from  elements  with  rela- 
tively high  mass  numbers.  However,  with  the  proper  selection 
of  lighter  materials  for  the  fuel  alloy,  container  and  reflec- 
tor, the  moderation  effects  May  Increase  to  the  extent  that 
this  fast  reactor  would  evolve  into  an  intermediate  reactor. 
It  is  believed  that  any  trend  toward  a  therml  reactor  would 
decrease  the  fuel  requirements  for  the  configuration. 

Conclusion 

Prom  a  nuclear  aspect  this  heat  exchanger,  which  was  com- 
posed by  volume  of  12.5  p«r  cent  alloyed  fuel,   12.5  per   cent 
molybdenum  structure,   and  75«0  per  cent  air,   was  critical  at 
a  diameter  and  length  of   kO  and  60  Inches,  respectively.     The 
total  weight  was  8,052  lbs  which  Included  3,622  lbs  of  enriched 


3U 

90  per  cent  uranium*  Approximately  $0   million  dollars  would 
be  required  for  thia  amount  of  uranium. 

The  flux  that  was  necessary  to  generate  a  maximum  wall 
temperature  of  2200*F  with  a  fuel  element  thickness,  0,0315 
in*  was  reasonable  at  10  ^n  per  sq  cm  sec* 

With  cost  as  a  major  consideration  there  is  probably  no 
point  in  a  fast  reactor  more  than  several  scattering  or  trans- 
port mean  free  paths  from  a  boundary.  An  intermediate  neutron 
energy  reactor  may  offer  the  possibility  of  reducing  the  fuel 
quantity  and  weight  for  this  configuration. 
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For  an  air  inlet  velocity  of  M  ■  0.2,  a  weighted  radial 
average  of  2000*1  'or  the  maximum  cylinder  temperature,  a 
frontal  area  of  7.7  sq  ft,  a  1*2  in  diameter,  and  a  60  In 
length,  the  amount  of  heat  transferred  was  substantial  enough 
to  produce  9,50"  lbs  of  thrust  with  a  flight  Hach  number  of 
2.$   at  37,500  ft. 

In  order  to  obtain  an  average  temperature  of  2000°l'  the 
fissionable  material  within  the  thin  right  oircular  cylinders 
had  to  reach  a  critical  condition.  This  required  3.622  lbs 
of  enriched  90  per  cent  uranium  with  an  approximate  cost  of 
50  million  dollars.   The  heterogeneous  nuclear  heat  exchanger 
was  critical  with  a  diameter  and  length  of  I4.0  and  60  in, 
respectively  and  a  total  weight  of  8,052  lbs. 

The  low  thrust  to  weight  ratio  may  be  compensated  for  by 
the  negligible  fuel  burn  up  when  compared  to  a  turbojet  sys- 
tem* 

Molybdenum  alloy,  0.1|5  per  cent  titanium,  which  was  pro- 
tected with  a  disilicide  coating  appeared  to  have  the  high 
temperature  properties  for  the  container  material.   Sandwiched 
between  the  molybdenum  was  an  alloy  of  uraniuwsirconium  which 
contained  79.3  per  cent  uranium  by  weight.   This  binary  system 
was  selected  because  of  the  2U50*F  melting  point  at  the  above 
composition.   Whether  a  metallurgical  bond  could  be  formed  be- 
tween the  molybdenum  and  fuel  alloy  was  undetermined. 


From  the  above  theoretical  conclusion*  and  limitin^. 
conditions,  the  over-all  feasibility  of  a  compact  axial  flow 
nuclear  heat  exchanger  was  favorable  for  thrust,  metal  avail- 
ability, and  nuclear  criticality,  borderline  for  thrust  to 
weight  ratio,  and  unfavorable  for  the  cost  of  fuel.  However, 
the  proposed  central  control  cor©  and  peripheral  reflector 
should  increase  thrust  and  reduce  the  cost  of  fuel  because  of 
savings  in  critical  mass. 

An  intermediate  neutron  energy  reactor  may  offer  the  pos- 
sibility of  reducing  the  fuel  requirements •  With  the  proper 
selection  of  lighter  materials  for  the  fuel  alloy,  container, 
and  reflector,  the  moderation  effects  s»y  increase  to  the  ex- 
tent that  this  fast  reactor  would  evolve  into  an  intermediate 
reactor. 
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EACH    ANNULUS: 
HYDRAULIC    DIA    =  D    -D, 
FLUID    FLOW    AREA  =  tt/4  (D22-D,2] 
FUEL    ELEMENTAREA   = 
rrt/2  (D,  +D2) 

HEAT   TRANSFER     AREA  = 
■n^D,*  D2)dx 


Figure  2.      Nuclear  heat  exchanger. 
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6000 


INLET   CONDITIONS: 
T0i2  =  884°  R 
pj«=  465  LB/IN2 


.3  .4  .5 

MACH  NO.   (INLET,  STATION    2) 
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Figure  3.  Thermal  choking  boundaries  for  a  constant  area 
frietionless  chamber  being  heated. 


Uo 


T0,AW        T0 


VT 


T0,AW~T 


x 


T 


0,W 


V/////////////y//////////////////////?//777^ 


he 
hm 


NO  HEAT  FLUX 
(ADIABATIC) 


HEATING 


Figure  k*     Relationship  of  various  temperatures  and  heat 
transfer  coefficients  where  TQ  -  T  *  V2/2gJC 
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Figure  6»     The  Brayton  cycle  for  temperature  and  entropy 

conditions  for  stations  corresponding  to  Figure  1. 
Values  for  dry  air  are  from  Reference   (7). 
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APPEHDIX  A 

Stepwise  Integration  for  Stagnation  Temperature 

and  Mach  dumber  as  a  Function  of  L/p  Ratio 

o 

The  procedure  used  for  integration  of  Equations  13  and 
H4,  are  outlined  in  Reference  (2)  and  arc  reproduced  in  part 
using  the  following  approximations: 

1.  The  integral  of  ATq/Tq   may  be  approximated  by 
2(T0.2  "  T0.1>/<T0,2  *  T0#1). 

2.  The  coefficients  of  dx/D  in  Equations  13  and  ll^  are 
taken  at  Mj  2  ■  (^   ♦  M2)/2. 

The  inlet  conditions  are  as  follows! 

»pR  •  0,6614. 

A   ■  1.32 
f   «  0.008 

^or  a  section  2  at  a  distance  x2  -  x1  downstream  from  section 
1,  the  following  computations  are  carried  out. 

a.  Let  <x2  -  ^i^e  ■  30 

b.  Estimate  Mp  ■  °»-*t  **  p  ■  °»2^ 

c.  Calculate  2(Ta  ,  •  ?A  .  )/(Tri  0  ♦  t   , 

0,2  0#1  '      <j#2  0fl)   from  Equation 

13. 


*! 

»  0.2 

T0,l 

*  88U*R 

po.i 

«  U6.5  lb/in2 

T0,w 

■  2U60»R 

•    *    io'« 


kl 


(To  r'oi1 

T0,2  ♦  T0#l 


O.w 


<*0.2  *X,iW 


-   1    + 


1  2        nl#2 


^       ii       ii  hi*      i 

1  \2  ♦  881; 


,»       2Wm.»  -  i  +    (1  -   .372H.2)(0fc>fl 
(T0,2        yW/2  1  ♦   .2  (.0625) 


2(1.32) (.008) (30) 


T0,2  "  l6lv0 


d.      Calculate  VU  from  Equation  15  and  compare  with  a»- 


>d  value  of  K 


M22  *   Ml2 


,T-^25),062?'     2(.008)(30). 


£< 


(1   ♦  l.k  x   .0625) (1 


.32)  [  .95] 


♦  2  x  I.I4.  x 


.0625 ( 


M 


!g     -  R^  *  0.0500 

M  2  *  0.0500  ♦  0.0U00 


l2 
**2 


*  0.3         l^ils  agrees  with  the  assumed  value. 


I 


1 


W 

e.  IT  atep  d  was  not  successful,  repeat  the  previous 
steps  with  an  improved  value  or  K2» 

f .  After  convergence  is  obtained,  the  process  ia  re- 
peated to  find  the  Wsch  nursber  at  a  section  3« 
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APPFUDIX  B 


Stepwise  Integration  for  Stagnation  Temperature 

and  t./D     Ratio  as  a  Function  of  Maeh  Number 

o 


The  procedures  used  for  integration  of  Equations  ll*.  and 
17  are  outlined  in  Reference  (2).  Equation  11*.  may  be  inte- 
grated approximately  over  a  short  Interval  of  Mach  number  by 
assuming  that  the  influence  coefficients  have  a  constant  value 
corresponding  to  (M.  •  M2)/2,  that  the  magnitude  of  TQ   is 
constant  during  integration  with  a  value  of  (TQ  ,  ♦  ffl  ?)/2t 
and  that  the  friction  factor  is  constant.   With  these  approxi- 
mations and  by  the  substitution  of  Equation  17  into  II4.,  Equa- 
tion Ik   integrates  Into 


n2     m  Hl 


2  {^ 
f0,l 


-  1) 


ft 


ft 


,1 


-f * 

T0,l     lQ9l 


^)-l 


(17a) 


let 


■  'or  the  same  Inlet  conditions  as  were  given  in  Appendix  A, 


M2  -  0,21+ 
Rj^  »  0.20 
\  z*  0922 


Kg 


0.0$76 
0,014-00 
0.0176 


v 


1 


So 

yrom  Reference  (2),  F-  ■  0.0^83  and  IL  ■  O.OO3I4B.   Equation 
17a  converged  after  three  attempts  to 


T 

J^£  «  1.31&  op  Tn  ?  »  1186*R 
x0,l  u»* 


Equation  1U  gave 


*2  "  *1  _  1     r„        f0.1  JB         1  Tn  2*785  -  1  _  1R  ft< 

"         fA    *    **  o      2T7o5oT  ^  2.785  -  1.3Wi      x8*85 


D.     2? 


Q.w    0.2 
^,1    *0,1 


Since  both  Mach  number  and  stagnation  temperatures  are  known 
at  each  section,  the  following  equations  were  used  for  comput- 
ing the  remaining  properties  of  the  air: 


T-   TQ   (1  ♦  dLSl  M.2) 
Tl    T0fl  (1  ♦  -Mp  K^) 


P2   M-   T 

p«|^  |^  ■  0,965 


V2   M    T 
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Pi    rl 

•• 

*2 

*0.2  m  ^2 
P0,l        h 

1   +  JLli  m  2 

1          2          2 

1  *  "^  Hl2 

■                                 « 

0.720 
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Tables  1-3  which  follow  have  calculate  values  Tor  the  various 
properties  of  air  as  a  i>unction  of  h/D     ratios  and  different 
oonstant  wall  temperatures. 


Table  1.   Air  properties  as  a  function  of  Kaoh  number  for  a 
constant  wall  temperature  of  266o0R  and  a  friction 
factor  of  •008 


H 

To 

r'/°e 

T 

P 

V 

e 

P0 

•R 

•R 

lb/in2 

ft/sec 

lb/ft3 

lb/In2 

0.20 

881* 

0 

880 

1*5.1 

291 

0.1385 

1+6.5 

0.21* 

1192 

11.90 

1182 

1*3.6 

1*05 

0.1000 

l*5.o 

0.28 

11+81* 

25.50 

11*66 

1*1.5 

$2k 

0.0770 

l*3.o 

0.32 

1738 

1*0.72 

1709 

39.2 

61*6 

0.0625 

1*0.6 

0.36 

1939 

56.01* 

1892 

36.8 

766 

0.0528 

38.1 

0.1*0 

2083 

69.91* 

2020 

31*.3 

881 

It  01*61 

35.6 

o.i*l* 

2182 

81.71* 

2101* 

31.8 

990 

0.01*10 

33.0 

0.1*8 

2252 

91.55 

2158 

29.6 

1097 

0.0372 

31.8 

0.^2 

2299 

99.19 

2182 

27.3 

1188 

0.031*0 

30.1* 

0,60 

2356 

109.52 

2182 

23.6 

1370 

0.029U 

28.6 

r 
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Table  2.   Air  properties  as  a  function  of  Maeh  number  for  a 
constant  wall  temperature  of  2l*6O0R  and  a  friction 
factor  of  .008 


M 

To 

L/D9 

T 

P 
lb/in2 

V 

ft/sec 

lb/ft3 

F0 
lb/in2 

0.20 

861v 

0 

880 

li.5.1 

291 

0.1385 

1*6.5 

0.21*. 

II8J4. 

13. U 

1180 

U3.U 

ttf 

0.0998 

UU.9 

0.28 

11*73 

29.32 

11*62 

bM 

527 

0.0769 

i*2.9 

0.32 

1711 

i*6.70 

1692 

38.9 

61*9 

0.0625 

2+0.3 

0.36 

1891 

63.98 

1863 

36.1* 

7*0 

0.0533 

37.6 

o,Uo 

2018 

79.33 

1970 

33*7 

878 

0.01*65 

3l*.8 

o.kk 

2100 

92.2S 

201*0 

31.2 

985 

o.oia5 

32.8 

0.1*8 

2155 

102.6 

2080 

23.9 

1085 

0.0376 

31.1+ 

0.52 

2192 

110.7 

2097 

26.8 

1182 

0.031*6 

30.2 

0.60 

2236 

121.9 

2100 

23.2 

1363 

0.0298 

27.8 

0.70 

2258 

128.2 

2050 

19.9 

1590 

0.0259 

25.8 

I 
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Table  3-   Air  properties  as  a  inaction  of  Maefa  number  for  a 
constant  wall  temperature  of  2060°R  and  a  friction 
factor  of  .008 


■ 

*0 

L/D# 

T 

P 

?0 

•R 

•R 

lb/in2 

lb/in2 

•20 

881+ 

0 

880 

1+5.1 

1+6.5 

•  2k 

1178 

17.95 

1160 

143.1+ 

1*4*6 

.28 

11*37 

39.75 

11+11 

i+1.0 

1+2.2 

.32 

1631 

63.55 

1592 

38.2 

39.1+ 

.36 

1760 

85.95 

1711 

IS.z 

36.3 

.1+0 

1839 

10U.7 

1775 

32.3 

33.3 

.10* 

1886 

102.0 

1811 

29.6 

31.1 

.1*8 

19114 

131.9 

1830 

26.9 

28.8 

.52 

1935 

11+0.6 

1830 

21+.3 

26.9 

•60 

191+5 

11+5.6 

1800 

21.1 

25.3 

' 
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APPFHTIX  C 

Calculation  for  Outlet  Stagnation  Temperatures  as  a 
Function  of  L/D  Ratioa  and  Maximum  Wall  Temperatures 

With  reference  to  Figure  2f  for  a  heat  balance  on  the 
fuel  cylinders  without  axial  conduction  requires  that 

Q    t 

To,w-To"-^r-    "^      •  <A> 

The  longitudinal  volumetric  heat  source  term  is  aero  when 
x  «*  0f  or  I*,  which  assumes  sero  extrapolation  distance. 

neglecting  axial  conduction  within  the  fluid  and  any  heat 
losses,  a  heat  balance  on  the  air  in  the  annull  requires  that 

f*CoD.  dT0 


To.w  "  *o  '  SlP  <ST    ■  <B> 


Equating  Equations  A  and  B  and  performing  the  required 
integration  results  in 


•o-'o.a-S^rfuj)  •[  i-"-P]- 


(c) 


Adding  Equations  A  and  C  to  eliminate  the  stagnation  tom« 
perature  of  the  air  yields 


4  * 


I*  •*0t*i*q 
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Saax     «*_  xx 


10,w   *0,2 


(bJ)[i-.o^p]     . 


• r^c *  i  *■ i  ■ *  -  cos~^  ■  *  (!° 


Substitution  of  Fquatlon  7  with  A  ■  1  into  Equation  D  results 
in 


t 


0,w   x0,2    C  C_ 


0  © 


(I) 


From  Equation  E  there  appears  to  be  a  maximum  which  repre- 
sents the  highest  temperature  of  the  fuel  cylinder  in  the 
longitudinal  direction*  The  point  at  which  the  maximum  is  at- 
tained is  found  by  differentiating  Equation  1  and  setting  the 
result  equal  to  zero. 


L 

*«.„~  *   r    are  tan 

max      % 


\~  h  ir )  <p) 


The  maximum  wall  temperature  occurring  in  the  fuel  cylin- 
der can  now  be  obtained  by  inserting  Equation  ■  into  equation 
E.  Yor   purposes  of  computation,  it  is  convenient  to  represent 
xx/L  by  the  symbol  a. 

a  S  xx/L  (G) 


( 
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By  substituting  G  Into  Equation  W9   results  in 

tan  a *  ~  ff  7T  (H) 


Utilizing  Equation  H,  Equation  E  reduces  to 


■"«x"'»-="^"  (t  )   I  '  -""«. 


(k)[ 


(I) 


From  Equation  17,  which  was  previously  derived, 

dT0     _  Pf  dx 

*0,w  -  10  u* 

and  introducing  Equation  A,  which  was  simplified  by  substitut- 
ing h  «  fGC  /2,  the  following  differential  equation  was  de- 
rived: 


2Q    t 

dTo-crfV    8l*-r  dx  <J) 

p  e 


Integrating  Equation  J  results  in 

US, 


To,3   To,2   gcx   V  r  J 


(K) 
P 

However,  Equation  K  May  be  combined  with  Equation  I  to  elirai* 
nate  (imAX   t,  which  results  in 


T 
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TQ,3  "  fM  *  2 


*  T3aX       * 

1  -  seca 


max 


(18) 


The   magnitude  of  Tn      utilized  the  previously  se- 


max 


lected  wall  temperatures.   Since  TQ  ~  was  *ke  given  inlet 

conditions,  a»  ,.  was  computed  for  various  I^L  ratios  result- 

ing  in  ffl  *  as  a  function  of  T0      and  <*^  •   Table  1  of 

*  max      ■■* 

Appendix  C  presents  calculated  values  of  TQ  ^  which  are  shown 

in  r'igure  $. 


Table  l*..   Outlet  stagnation  temperatures  as  a  function  l*/D 
ratios  and  maximum  wall  temperatures 


L/D 


l|-0 


60 


80 


100 


120 


11*0 


tana 


max 


seca 


1-seca 


'max 


f 


o»3 


o#3 


-i>9lo  -3*270  -2.1£o  -1.965  -1.635  -i.l*o5 

-5.010  -3.1420  -2.650  -2.210  -1.920  -1.728 

6.010       fc*kH       3.650       3.210       2.920       2.726 

TQ  2  «  88ir»R  and  ?Q             ■  266 0eR 
* *    max 

•R       1U75         1688         1858         1990         2101         2181*. 

TQ  2  *  88U#R  and  TQ  w         «  2l*©0»R 

*  *   max 

•R       li*08         1597         1714B         1865         1961*.         2038 

T0  2  *  88if*R  and  TQ  a  206 0«R 

9  'max 


•tJ 


•R       1275         XkV. 


1529         16X6         1689         l7l|5 
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APPE*I>IX  D 
Critical  Calculations  for  a  Homogeneous  Unreflected  React or 


The  effective  neutron  en&r&j  was  taken  as  0.3  Hev  based 
on  flux  spectra  computed  for  existing  fast  reactors  by  Murray 
(18#  p.69).  From  the  same  Reference  the  following  constants 
were  used: 

XT^   8  <T&   «  0,18  barns    0"f  »  0.0         0"t  *  7.1+  barns 

tr-*^  j  <T '    m   1.78  barns    0\.  «  1.5  »  barns    T^  m   7.I4.  barns 

V»  2.51 

From  Hughes  and  Harvey  (19)  at  0.3  Mev,  the  total  cross 
sections  of  the  selected  reactor  materials  and  air  were  as 
follows : 

Zr  :  0"T  *  9.0  barns 

Mo  t  <TL  »  9*0  barns 

Jf  1  <T  *  l+.o  barns 

0  1  0"T  «  3.0  barns 

The  above  total  cross  sections  were  considered  to  be  equal  ap- 
proximately to  the  scattering  or  transport  cross  sections. 

A  90  per  eent  enriched  uranium  mixture  by  weight  produced 


235  »k 


238  N8  ♦  Z3S  ns 


SB      .90 


or 


*a  m  0.105 


otf*    «ff 


The  solution  to  Equation  22  with  the  above  result  yielded 

*7  ■    7= r^*| *  2.51  r~7r — ^A^w,^  *  2'°9  • 

'      <Ta,5   *V,8  (Hq/m^)      1*?G  *  M8)U05) 

Since  the  alloyed  density  of  the  uranium- zirconium  system 
was  unknown,  the  alloy,  79  »3  per  cent  uranium  by  weight,  was 
treated  as  a  mixture  to  obtain  the  following  volume  fractions: 

f238  m   °*°5* 
f23«;  3  0.512 

f2r  "  0mk3Z 
A   hydraulic  diameter,  D  »  0.5  in,  was  obtained  with  a 
selected  r,/r  »  120  and  L  •  60  In.   The  following  thicknesses 
were  assumed! 

Uranium- zirconium  alloy:  0.  "53125  in 

Double  wall  of  molybdenum:  0.03125  in 

Hydraulic  diameter/2:  0.25000  in 

Thickness  of  air  space:  0.18750  in 

The  approximate  component  volume  fractions  obtained  from  the 
above  solidity  factors  were  as  follows! 


r!f 

SB 

0.5925 

*0 

St 

0.1575 

f23Q 

at 

0.0070 

f235 

3S 

0.06I4.0 

Zr 

at 

0.0514.0 

fK0 

at 

0.1250 
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The  macroscopic  cross  sections  were  computed  Tor  each  of 
the  above  volwse  fractions  which  gave  for  a  total 


£.  «  0.1163  era*1    and    2   *  0.0056  cm"1 


Solution  to  "mat  Ion  27   yielded 

B2  m   21.7  x  10"^  era"2   • 

With  a  physical   length  equal  to  60  in  or  tfteJt  era  and  the 
transport  mean  free  path  equal  to  8.U5  cm,   the  length  includ- 
ing the  extrapolation  distance  was 

H  •  h  +   .7lAt  «  152.14-  •  6.0  ■  156.U  o« 

The  solution  to  Equation  29  for  the  critical  radius  in- 
cluding the  extrapolation  distance  was  56.7  em;  the  physical 
radius  was  50.7  cm  or  20  in. 

Since  the  total  volume  was  12.35  x  Mr    eu  cm,  the  follow- 
ing component  weights  were  determined  by  using  the  previously 
computed  volume  fractions  with  corresponding  densities; 


U*38 

362  lbs 

U235 

3260  lbs 

ftp 

960  lbs 

Mo 

3U70  lbs 

,  Weight 

8052  lbs 

■iro-r. 


i   i^M 


aiJ 


m  uv? 


?v  $a»m 
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Equation  Kf   Appendix  C,   the  flux  that  was  necessary 
to  generate  a  Bjaximum  wall  temperature  of  22     °     with  a  fuel 
element  thickness,    0«03l£  in,   was  obtained  with 


Smi 


Z|£  £f  t  (r-/oe) 


where 


E  -  200  Mev/f  lesion 

Mev  ■  1#$2  x  10"16  BTU 

G  *  kO.k   lb/ft2  sec 

C  ■  .2671  BTU/lb  • 
P 


t/T>  *  120 

© 

T0  3  "  T0  2  *  2101  -  88^*  frOW  P1£ure  *• 

The  flux  obtained  frosi  the  above  equation  was  8#25  x  10  Tn/ 


sec. 
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